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Abstract 

J>,- We study the distribution of lepton pairs from the second hghtest 

D I neutrahno decay X2 ~^ Xi^^^- This decay mode is important to mea- 

"^ ' sure the mass difference between X2 ^^'^ ^^^ hghtest neutrahno Xii 

which helps to determine the parameters of the minimal supersym- 
k> , metric standard model at the CERN LHC. We found that the decay 

H I distribution strongly depends on the values of underlying MSSM pa- 

rameters. For some extreme cases, the amplitude near the end point 
of the lepton invariant mass distribution can be suppressed so strongly 
that one needs the information of the whole niu distribution to extract 
m^o — m^o. On the other hand, if systematic errors on the acceptance 
can be controlled, this distribution can be used to constrain slepton 
masses and the ZX2X1 coupling. Measurements of the velocity distri- 
bution of X2 from samples near the end point of the mu distribution, 
and of the asymmetry of the pT of leptons, would be useful to reduce 
the systematic errors. 



1 Introduction 

If supersymmetry is realized in nature, it promises exciting possibilities for 
future collider physics — the discovery of sparticles. If the scale of super- 
symmetry breaking is around 1 TeV (as preferred by fine-tuning arguments 
related to problems in the Higgs sector), many sparticles will be produced 
at future colliders such as the Fermilab Tevatron upgrade, the CERN Large 
Hadron Collider (LHC), or future e^e~ colliders proposed by DESY, KEK, 
and SLAC. Sparticles produce unique signatures in the detectors and will be 
discovered quite easily. 

It has been pointed out that one cannot only discover those sparticles, but 
can also study their detailed nature in future e'^e~ colliders. By measuring 
sparticle masses, production cross sections for a polarized electron beam, 
and other distributions, we will measure soft supersymmetry breaking mass 
parameters |]l], ||, H, ^, |] and prove supersymmetric relations 0. Studies 
at future e^e~ colliders should reveal details of the mechanism to break 
supersymmetry. 

Corresponding studies for the LHC have been performed in [^ Q in the 
framework of the minimal supergravity model. These analyses show that 
a precise determination of the model parameters is possible; the LHC is 
especially powerful when the gluino decay g ^ bb followed by 6 — *> bx2, X2 ~^ 
l~^l~Xi can be identified by tagging the bottom jets and the hard leptons. 
One of the key tricks of the studies is the measurement of the end point 
of the invariant mass distribution of the lepton pair with same flavor and 
opposite charges. The end point determines the mass difference between the 
second lightest and the lightest neutralino [^, rriy^o — m^o. The kinematically 
constrained nature of the end point samples allows one to reconstruct the 
decay chains to determine the parameters of minimal supergravity model 
U^, 1^. At the end, all or some of the parameters in the model would be 
determined.^ 

However, the minimal supergravity model is not the only model for spon- 
taneous supersymmetry breaking in a phenomenologically consistent manner. 
Moreover, there are various possibilities even within the supergravity model, 
as summarized in Ref. |]12[ for example. Other mechanisms to break super- 



symmetry (SUSY) are also discussed extensively |Tl], |13 . 

Assuming that the LHC can determine all parameters within the mini- 
mal supergravity model by looking into some signal distributions, the next 
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"'^A similar analysis can be do ne PQ I for the parameter dependence of gauge mediated 
supersymmetry breaking models 



question is what we should look into to over-constrain the model, or to de- 
termine deviations from the minimal supergravity model. As we mentioned 
before, studies in this direction have been done in great detail for future e^e~ 
colliders. For hadron colliders, such a study would be complicated because 
more than one SUSY channel can contribute to a particular class of events 
since different sparticles are simultaneously produced. 

In this paper, we discuss the parameter dependence of the distribution of 
the three body decay X2 ~~^ Xi^~^^~- The decay branching ratio depends on 
parameters in the minimal supersymmetric standard model (MSSM) rather 
sensitively due to (possibly negative) interference between slepton and Z° 
exchange contributions ||l^ |^, |T6|, |17|. We point out that not only the 
branching ratio but also the decay distribution is sensitive to the MSSM 
parameters, giving us an extra handle to determine mj, tan/3. Mi and fi, 
independent of the SUSY breaking mechanism. 

The organization of this paper is as follows. In section 2, we discuss the 
MSSM parameter dependence of the lepton invariant mass distribution aris- 
ing from X2 decay. Near the region where interference between Z^ exchange 
and / exchange becomes important, the decay distribution is very sensitive to 
rrii. We point out that studying the whole distribution of the lepton invari- 
ant mass mil can be very important even for extracting the end point; events 
near the end point could be so few that one could misidentify it without this 
information. We also discuss the chiral structure of the amplitude, and show 
an interesting parameter dependence of the tau polarization in X2 — »• t~^t~Xi 
decays. 

In section 3, we discuss effects of cuts on observed distributions. Though 
mil of lepton pairs from x^ decay is independent of the X2 boost, each lep- 
ton energy depends on the parent Xo, momentum. Therefore we expect two 
apparent sources of systematic error: (A) uncertainty of the lepton energy 
distribution in the X2 ^^^^ frame, and (B) uncertainty of the parent X2 ^^O" 
mentum distribution. It would be helpful to reduce these two errors in order 
to maximize the physical information that can be extracted from the ob- 
served m^ii distribution. We argue that a measurement of the asymmetry of 
lepton energies would reduce the systematic errors from (A), while the lepton 
energy distribution near the end point of the mu distribution constrains (B). 
In section 4 we study if the mu distribution constrains mj and the param- 
eters in the neutralino mass matrix. Section 5 is devoted to discussion and 
comments. 



2 Lepton Invariant Mass Distribution 

The branching ratio of the leptonic decay of the second hghtest neutrahno, 
X2 ~^ Xi^^^~; is known to be very sensitive to the values of the underlying 
parameters. Three body decays of X2 ^^^ dominant 0, |15|, |l6l as long 
as X2 — *■ Xi^°) X2 ~^ ^^ a^6 not open. The dependence is enhanced by 
the negative interference between the decay amplitude from Z^ exchange 
and that from slepton exchange. In the minimal supergravity model, the 
interference would be maximal when slepton masses are around 200 GeV 
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In this section we show that the effect of the interference appears not 
only in the branching ratios, but also in the decay distributions, such as the 
distribution of the invariant mass mu of two the leptons. 

We first show the differential partial width of the decay Xa ~~^ X%ff fo^^ ^ 
general light fermion /. We assume that the mass and Yukawa coupling of /, 
and left-right mixing of / are negligible. The decay amplitude then consists of 
the Z^ channel and the f^j^ channel. The squared, spin averaged^ amplitude 
|A^P of the decay Xa(p) ~^ XBiP)fi^)fil) is written as follows: 

\M\' = 2{Al, + Al^){l-y){y-rlJ+2{Al^ + Al,){l-x){x-rlJ 



-4:{AllArl + ARRALR)r^j,z, 



with 
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All - -9z 2 Z^ o5'2 

Arl = -All^V ^ x), 

A ^^,2'^BA'^R I Y^ -'-^2"AX"BX 
^LR — -Qz 2~ + Z^ o5'2 2~' 

2 z — ri, ^2 X — rj 
Arr = -ALRiy^x), (2) 

where .[') = T^jl - QfS^w. 4^ = -Qfs'w, with si, = .^H^. The forms 
of the x^x"^*^ couplings z^j^ and x"// couplings (o-ax; ^ax) ^^^ given in 
Appendix, (x, y, z) are phase space parameters of the decay defined as 

X = (p + g)7m|,^, y={p + qf/m\o^, z = {q + qf / m\o^, (3) 



^The polarization of X2 might affect observed decay distributions. The effect depend 
on the process for X2 production, and is not discussed in this paper. 



with x + y + z = l + rn^-o /m?-o . For convenience, we have introduced mass 
ratios r^^, rz, and r? as 

^XB = "^x°b/"^xV ^^ = rnz/m^o^, rj^ = mj^/m^o^. (4) 

Note that {All,Arl) involve {fhjR) while {Alr^Aur) involve (fRji)- 
They do not interfere for mj = 0. 

The partial decay width is given by 

^^ ',-,0 . ~0 rn-, _ ^C _ I . .,2/ ., ^ _ 1 , ^2 



where A^^ = 3(1) for / = q{l). The range of (x, y) is given by the conditions 

z{xy-rl ) > 0, 



r'^g <x < 1, 
rL < 1/ < 1, 
X + 1/ + 2; = 1 + r?^. (6) 

Now we consider the case of X2 decay into Xi^^^~ under the assumption 
that all two body decays of X2 ^^^ kinematically forbidden.^ 

In the phase space of the decay X2 ~^ Xi^^^ > the Z° exchange amplitude 
and / exchange amplitude behave differently. When the Z^ contribution 
dominates, distributions are enhanced in the region of large mfi = zm^-o. 

In contrast, when the / exchange contribution dominates, distributions are 
enhanced in regions with large x and/or large y, therefore in small mn and 
large |ii^[-!^* — El%^^\ region. Here 

Elf' = (1 - x)m^o/2, i^r* = (1 - y)m^o/2, (7) 

are lepton energies in the X2 ^^st frame. 

We consider the case where 2Mi ~ M2 ^ |/^|) a typical case in the 
minimal supergravity model. In this case, X2 is Wino-like and Xi is Bino- 
like. An interesting property in this case is that the Z^ and / amplitudes 
could be of comparable size in some region of phase space. Furthermore, 
their interference is generally destructive for leptonic decays. These effects 
cause complicated situations, which we discuss below. 

We show numerical results for the mu distribution. For illustration, we 
use two sets of parameters for the neutralino sector, (A) and (B), shown in 
Table ^ These values are fixed to give the same masses for three inos, {m^o , 

■^In some region of parameter space, we may study the three body decays of X2 even if 
two body decays are open p8| . Such a study is beyond the scope of this paper. 



set 


Ml 


M2 


/i 


tan/? 


(A) 


70 


140 


-300 


4 


(B) 


77.6 


165 


286 


4 



Table 1: Parameter sets for neutralinos. All entries with mass units are in GeV. 

m^o, m-+) = (71.4, 140.1, 320.6) GeV. For calculating the branching ratios, 
we take generation-independent slepton masses and a universal soft SUSY 
breaking squark mass mq = 500 GeV. 
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Figure 1: Invariant mass distribution of the lepton pairs from X2 three body decay. The 
neutralino parameters are taken as set (A) (see Table w), and universal slepton masses 
mj = mj^ = mj^ are 170 GeV (thick sohd), 220 GeV (dashed), 270 GeV (dotted), 320 
GeV, and 500 G^eV. For m^ = 500 GeV, Br{x2 ^ e+e-^?) =11%, 9.5%, 4.1%, 2.1%, 
1.9%, respectively. The total number of events of each curve is lO'*. 

In Fig. 1, we show the mu distribution of the decay X2 ~^ ^^^Xi for 
parameter set (A) and varying mj from 170 GeV to 500 GeV. 

Because m^o and m^o are fixed, the end points of the distributions mfl^^ = 
68.7 GeV are same for each curve, while the shape of the distribution changes 
drastically with slepton mass. For a slepton mass of 170 GeV (thick solid 
line), the decay proceeds dominantly through slepton exchanges, therefore 

^Heavy sleptons and Bino-like Xi is cosmologically disfavored because it leads to a large 
relic mass density [|l9[. However, this constraint can easily be evaded if Xi can decay, or 
if there is late time entropy production. 



the mil distribution is suppressed near m™^^. On the other hand, once slep- 
ton exchange is suppressed by its mass, Z^ exchange dominates and the 
distribution peaks sharply near m™'^^. Notice that the leptonic branching 
ratio decreases as mj increases, but remains larger than ~2% for the values 
of mj taken in the figure. 
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Figure 2: mu distribution for parameter set (B) (/x > 0). 



arc 170 GeV (thick 



solid), 220 GeV (dashed), 250 GeV (dotted), 370 GeV, and 500 GeV. The branching ratio 
Br{x2 -^ e+e^xl) for m^ = 500 GeV is 6.6%, 2.9%, 0.9%, 1%, 1.8%, respectively. The 
thick dashed line is an example with complete cancellation of the amplitude near the upper 
end point. The branching ratio to e+e^Xi is 1.8% for this case. The total number of the 
events is 10"* for each curve. 



In Fig. 2, we show an example for /i > 0, parameter set (B). The de- 
pendence on the slepton mass is different from the previous case. As mj 
increases from 170 GeV, mu distribution becomes softer. For m^j > 250 GeV, 
a second peak appears due to strong cancellation of Z^ exchange and slepton 
exchange contributions for a certain value of m^u. At the same time, the 
branching ratio reaches its minimum at m^j ~ 300 GeV, much less than 1%. 
For m^i ^ 370 GeV, it increases again above 1%. 

Notably, one can find slepton masses where a complete cancellation occurs 
very close to the end point m™*^^ of the m,ii distribution. The thick dashed 
line shows distribution for mr = 245 GeV and ttit = 204 GeV. Events near 
the end point {m,f^^ — m^u < 4 GeV) becomes too few, and it is very hard to 
observe the real end point for this case. 

The lepton invariant mass distribution is an important tool for studying 



supersjTiinietric models at hadron colliders. In |p, a case study is done in a 
scenario where decays g ^ bb and b -^ bx2 have substantial branching ratios. 
X2 production is enhanced by the large gluino production cross section, and 
the S/N ratio could be improved substantially by requiring three or four 
bottom jets in the final states. For Br{x2 ~^ e"'"e~Xi) ~ 16 %, S/N goes 
well above 10. Even if gluino decay into bb is closed, the lepton invariant 
mass distribution can be measured (case 4 of p). In this case, one can 
subtract most backgrounds using lepton pair samples with opposite charge 
and different flavors. 

The most important aspect of the lepton invariant mass distribution in 
these studies is the determination of the end point, which is expected to coin- 
cide with m^o —m^o. Furthermore, choosing events near the end point helps 
one to relate the velocity of the two lepton system (3ii to that of the neutrali- 
nos (3^0 and /9^o, so that one can reconstruct the cascade decay chain. Errors 
on the mass difference of 2% ~ 0.1% are claimed depending on statistics p. 

The slepton mass dependence of X2 decay distributions suggests that not 
only the end point of the distributions but also the distributions themselves 
contain information about the underlying parameters such as mj. The neg- 
ative side of this is that the fltted end point may depend on the assumed 
values of these parameters, introducing additional systematic errors in the 
flt. For some extreme case shown in Fig. 2, the observed end point of the 
lepton invariant mass distribution does not coincide with m^o — m^o. Note 
that realistic simulations including the parameter dependence of the decay 
distribution are not available for hadron colliders so far. In the commonly 
used Monte Carlo (MC) simulators ISAJET ^ and SPYTHIA ^, the 
three body decay distribution of sparticles is approximated by the phase 
space distribution, while branching ratios are calculated by full expressions.^ 
See section 3 (especially Fig. 4) for comparison between real mu distribution 
and the distribution in phase space approximation. 

For the case shown in Figs. 1 and 2, Br{x2 -^ Xi^^^) is not always 
large. It has not been yet studied systematically if it is possible to mea- 
sure the lepton decay distribution at future hadron and e~^e~ colliders when 
Br{x2 -^ Xi^~^^~) is small. Given small Standard Model backgrounds, a 
study of distribution may be possible when Br{x2 ~^ e"'"e~Xi) ~ 2% for not 
too heavy gluino. At future lepton colhders, it is very easy to get a clean sig- 
nal. However, the production cross section is rather small for heavy sleptons. 



^The most recent ISAJET release (ISAJET 7.43) allows to simulate the effect of exact 
matrix elements for all three body decay distributions. In the codes used for CERN LEP 
[p2[ and JLCl [ESl studies, the effect is already included. 



At e~^e~ collider with large luminosity {J C > 500 fb^^/yr) [^, the region 
with cr(e"'"e~ -^ X2X2) > 20 fb may be studied with samples containing more 
than 1000 events, even if Br{x2 -^ e+e^x^) ~ ^ few %. 

In Figs. 1 and 2, we found that the decay distributions depend on sign(/i). 
The difference comes from neutralino mixings. For the parameter set (A), X2 
has small Bino component; (iVai, N22, N23, N24) = (7.25x 10-^ 0.957, -0.281, 
-6.90 X 10^2) „here x- = N^B + Ni2W^ + Ni^H^ + Ni^Hl It does not 
couple to Ir effectively [see Eq. (|A^) ], therefore decay into Ir proceeds dom- 
inantly through Z^ exchange. In Fig. 3(a) we show the niu distribution 
for Z"*"/]^ (dotted lines) and /"""/j^ (dashed lines) for m^- = 220 GeV(thick lines) 
and 270 GeV(thin lines). When mi increases from 220 GeV to 270 GeV, 



\M, 



\M.{xV'^^l)\'^ is suppressed near mfl^'^ and also becomes smaller 



due to interference. The relative importance of \M.r\ increases, leading to 
a total distribution that is more strongly peaked near m™^^. On the other 
hand, for parameter set (B) (yU > 0), X2 still has substantial Bino component; 
(iV2i,iV22, iV23,iV24) = (0.196,0.902,-0.322,0.212). As m^- increases, nega- 
tive interference appears near mu = mf^'^^ for both \Ml\ and |A^_r|, leading 
to a suppression of the mu distribution near mf^^'^. 
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Figure 3: (a) mn distribution from x^ 



Xil 



decays for mr — 220 GeV and mr = 270 



GeV and parameter set (A). The dashed hnes show Z+/£ distribution, the dotted hues for 
/+^^ distribution, the sohd hnes for totaL The curves are normahzed so that the total 
widths are the same, (b) Polarization of lepton l~ [= {Nr — NL)/{Nfi, + Nl)] arising 
from X2 decay. Average polarization is Pave = —0.839 for mj = 220 GeV and —0.428 for 



mj = 270 GeV. 



In Fig. 3(a), negative interference reduces the amplitude near m]^^^ for / 



L' 



while the amplitude for Ij^ is increased. This causes a strong mu dependence 
of the lepton polarization P;. In Fig. 3(b), we plot this polarization as a func- 
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tion of mil for different slepton masses. They differ dramatically near mfl^^. 
In neutralino decays into t^t^Xi^ Pt may be observed through the decay dis- 
tributions of the r leptons. In r"*^ ^ p"*^ ^ tt'^tt", ai —>■ 7r^7r*^7r° decays, the 
E^±/Ejet distribution depends on the parent r polarization drastically P5| . 
Experimentally, the fine momentum resolution [|I], |^ of the detector could 
be used to identify these tau decay products |^ at future e~^e~ colliders. At 
hadron colliders, the momentum of charged tracks and the energy deposited 
in the electromagnetic calorimeter would give the same information. Impli- 
cations of T polarization measurements in analyses of new physics have been 
discussed for charged Higgs bosons produced at the upgraded Tevatron and 
LHC colliders p6|, and for f at e~^e~ colhders ||2^, ||. 



Due to the missing tau neutrinos, one would not be able to measure the 
invariant mass of two tau leptons. Nevertheless the Pr dependence on the 
invariant mass of the two r jets might be seen in future collider experiments. 
Note that in r ^ p or ai decays, the final vector meson carries a substantial 
part of the parent r momentum, therefore the smearing of the distribution 
is less severe than for decays into vr^, /x, and e. 

Several comments are in order. In Eq. (|I|), we have neglected the Yukawa 
couplings of leptons and slepton left right mixing. For f, these effects could 
be very important if tan/? is large. Notice that their leading contribution 
fiips the chirality of the r lepton |2^ . For three body decays, studying the 
correlation of two tau decay distributions would reveal the helicity flipping 
and conserving contributions separately. 

In most numerical calculations in this paper, we assume universality of 
slepton masses. However, in supersymmetric model, staus could be lighter 
than the other sleptons for various reasons. The running of stau soft SUSY 
breaking masses from the Planck scale ||2^ and stau left-right mixing |29 



could enhance decays into rf or t'^t^Xv Experimental consequences of 
such scenarios have recently been widely discussed fS^, 0. Also models 
with lighter third generation sparticles have been constructed to naturally 
avoid the flavor changing neutral current problem [jl3|. The two body decay 



branching ratio into rf or the three body decay branching ratio into ttXi 
and the decay distribution might be different from those for leptons in the 
flrst two generations. The study of the X2 ~^ t'^'t'Xi decay in addition to the 
other leptonic modes could be an important handle to identify such models. 



3 Correcting Acceptance Errors 

In collider experiments, we need cuts to reduce backgrounds, and they sub- 
stantially change distributions we are interested in. This effect can be cor- 
rected by Monte Carlo simulations once model parameters are fixed. But it 
is still worthwhile to investigate how the distributions are modified by the 
cuts, and if this effect can be estimated in a model-independent way. 

Let us consider X2 decays at the LHC, where X2 comes from the cascade 
decay of heavier squarks and gluinos: 

f gq (if rrig > rrig), 
^ I QX'2. 

~g - j ^^:*.°1 . (8) 

[ qq (if rriq < m-g), 

and X2 ^^^y decay further into l^l^x\- These decay processes have been 
shown to have small backgrounds. In j^, ^ decay into hh followed hyh ^ 6x25 
and g decay into qqxH were studied. After lepton transverse momentum cuts, 
cuts on total transverse energy Et and total missing transverse energy i^y, h 
tagging (for the former), and subtraction of backgrounds estimated from the 
opposite charge - different flavor lepton sample, the irtu distribution from X2 
decays can be measured over a wide region of mu . 

The observed distributions are modified by the cuts, and they depend on 
the decay processes through which x^ is produced. One might worry that 
those cuts are strongly correlated with mu, i.e., particular regions oiirtu have 
large (small) acceptance, therefore the whole distribution may become quite 
different from the original one. Indeed, the cuts on lepton energy should 
affect the mu distribution directly. On the other hand, cuts on the total Et 
and EJt should be less correlated with the momenta of the leptons. In the 
following, we therefore discuss the effect of cuts on lepton energies on the 
lepton invariant mass distribution. 

In Ref. [§, a lepton transverse momentum cut, plj. > 10 GeV ~15 GeV, 
is applied to all leptons. The m,u distribution after this cut should depend 
on the El distributions of X2 decays in the X2 ^^^^ frame, and also on the 
X2 momentum distribution through the boost of leptons. In Fig. 4, we plot 
the m^u distribution for our standard parameter set (A) and m^- = 250 GeV, 
requiring that the lepton energy Ei in the X2 ^^st frame is larger than (thick 
solid), 10 GeV (thick dashed), 15 GeV (thick dotted), respectively. We also 
plot the corresponding distributions in phase space approximation. 

Figure 4 reproduces the quantitative effects of pip cut. Because of the 
lepton pt cuts, leptons with Ei < p^"* have no chance of being accepted 
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Figure 4: Invariant mass distributions of the Icpton pairs from X2 decay for the parameter 



set (A) and 



250 GeV (thick hncs). Ei > (solid hne), 10 GcV (dashed), and 15 



GeV (dotted) is required. We also show the distributions in phase space approximation 
by thin lines. 



unless the parent neutralino has enough transverse momentum. Therefore 
the dashed and dotted hues show quantitative nature of the mu distribution 
after the pip cut when the momentum of the parent neutrahno is neghgible. 
Typically, events with large mu have a better chance of being accepted. 
We can see that both of the leptons have energy larger than 15 GeV if 
mil > 54 GeV. Ei^^^ ~ 35 GeV in the region close to m™'^^. Unless leptons 
go in the beam direction, they would be accepted. Although the actual 
distribution after the pip cut would depend on the px distribution of X2! it is 
still qualitatively true that events with small m^u would be more affected by 
the cuts. Note that, whenever a large cancellation of the amplitude occurs 
at the mil end point, the m^u distribution near the end point also differs 
substantially from those without cancellation. The acceptance near the end 
point should be large and should depend on m^u only weakly, so that we could 
distinguish a "fake" end point from the kinematical one through the study 
of the distribution. 

In Fig. 4, the mu distribution is softer in the phase space approximation. 
However, when Ef*^^^ > 10 GeV is required, the number of events with smaller 
m,u is reduced significantly compared with that for parameter set (A). The 
qualitative difference between the two curves becomes less significant after 
■rest ^ ]^Q QgY jg ];equired. 



rarest 
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Figure 5: Ai = |£;[f * - ^[5''*|/ (£;[f' + S™''*) distribution of the lepton pair from X2 
decays with 15 GeV< mu < 40 GeV and Ei > 10 GeV. The thick sohd hne is for nij = 250 
GeV and the parameter set (A), the dashed hne for mj =170 GeV and the parameter set 
(B), and the sohd hne for phase space approximation. Gurves are normahzed to coincide 
at Ai = 0. 



The difference of the acceptance in smaller mu region can be seen when we 
plot lepton energy difference. In Fig. 5 we plot Ai = \Ei+ — Ei- \/{Ei+ + Ei~) 
in the rest frame of X2 requiring Ei > 10 GeV and^ 15 GeV < mu < 40 
GeV. The distribution reaches its maximum a.t Ai = for the parameter set 
(A) with m,j = 250 GeV (thick line), while for phase space approximation 
(solid) the Ai distributions are roughly fiat. This is because the amplitude 
near Ei'^^^ = is suppressed, and agrees with the quantitative difference of 
the acceptance in small mu region found in Fig. 4. The acceptance should 
also depend on underlying MSSM parameters. We show the distribution for 
the parameter set (B) and mj =170 GeV, which has similar m,u distribution 
to yU < case. Due to the small slepton masses, the amplitude is enhanced 
in smaller Ei region, therefore the Ai distribution is more fiat than for the 



m^i = 250 GeV case. This suggests that Ai and acceptance may depend on 
slepton masses when X2 ~^ ^^ is about to open. 

It might be possible to determine the lepton energy asymmetry directly. 
In Ref. [01, the correlation between the Ai and Al[ distributions is studied 



^Hcrc we restrict mu from 15 GeV to 40 GeV. We may have a background of 7* -^ l^l^ 
in small mu region. For tuu near the end point, leptons tend to have equal energy, and 
the contribution to Ai is insignificant. 
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for the X2 decay into IrI. Here Af = |p^ — p?^ |/(Pt + Pt ); where pip and 
pip are transverse momenta of the two identified leptons. In the two body 
decay, Ai in the X2 ^^st frame is restricted to a rather small region. Although 
X2 originates from gluino or squark decays and should have some transverse 
momentum for the case studied in that paper, the Af distribution is still 
peaked where the Ai distribution is located. This shows that Af reflects the 
decay distribution in the X2 ^^^^ frame. 

The acceptance would also depend on the X2 momentum distribution. 
If X2 h^s large transverse momentum, events with smaller Ej^^^ would have 
a better chance of being accepted. Notice that if i?[^^* < p^"*, this is the 
only possibility for the events being accepted; certainly the X2 transverse 
momentum distribution must be determined somehow to make a reliable 
estimate of the acceptance. 

Naturally, an average pj.^o of the order of 'mg(^g) is expected when the X2 
comes from q or g decay. The tree level kinematical distribution is deter- 
mined by parent and daughter masses if the X2 originated from a two body 
decay. Existing MC simulations describe the two body decay distributions 
correctly. For the three body decay, the distributions are different from those 
in phase space approximation. However, there would be no significant pa- 
rameter dependence because only squark exchange contributes to the decay 
modes in Eq. (j^). If we have direct information on those masses, it may be 
rather easy to estimate the transverse momentum distribution^ of xt- 

The X2 distribution can be measured directly if there is substantial statis- 
tics near the mu end point. At the end point, the total momentum of the 
lepton pair is zero in the X2 ^^st frame, and we could reconstruct the velocity 
of X25 Px°y from the momentum of the lepton pair [0, |]. This X2 distribu- 
tion may be convoluted with the mu distribution and the Ai distribution in 
the X2 '^^st frame to obtain the corresponding distributions in the laboratory 
frame. Even though the observed decay distributions are sensitive to the p^^o 
spectrum, they might thus be corrected by measuring the (3^0 distribution 
using events with mu near its end point. 

So far we have only discussed the invariant mass distribution of lepton 
pairs from X2 decays, which is known to be important for SUSY studies at 
the LHC. In contrast, for the Tevatron upgrade, the trilepton mode is more 
important to discover supersymmetry [^, |30|. This signal comes from the co- 
production of X2xf s-iid their decays into leptons. For this particular mode, 
the energy distribution of leptons from chargino decay must be considered 

^However, the transverse momentum of gluino receives (SUSY) QCD corrections. Cor- 
rect estimation of the distribution would be very important. 
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in addition to the neutralino decay distributions. 
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Figure 6: Energy distribution of the charged lepton from x^ — > Xil^'^ decay in the xt 
rest frame. We take parameter set (A) (/x < 0) and (B) (/i > 0) and vary mj as indicated 
in the figure. The thick sohd hue shows the distribution in the phase space approximation. 
The leptonic branching ratio of xt is above 8% for the parameters used in the figure. 

In Fig. 6, we show the lepton energy distribution of chargino xt ^^ the 
chargino rest frame for the parameter sets (A) and (B). The amphtude of 
Xi^ decay is easily obtained by replacing relevant masses and couplings in 
Eq. (|l]) . The parameter dependence is not very strong if the slepton mass is 
much above 200 GeV, but if it is close to xt: the lepton energy distribution 
sensitively depends on the slepton mass. The branching ratio to the lepton 
is also larger in this region making the trilepton mode more promising. For 
parameter set (A), the distribution is almost flat in Ei for rrij = 170 GeV. 
Such a dependence of the energy distribution in the xt rest frame on MSSM 
parameters will affect the observed distribution of leptons in the lab frame. 
Notice that we must also pay attention to the effect of other cuts, e.g., on 
the total Et or ]^t- For example, for X2xf co-production, the two x have 
balanced transverse momentum unlike in the case where X2 comes from g 
decay. The two Xi from X2 and xf decays tend to be back to back to 
each other. This makes the total I^t smaller, especially when the lepton 
energies are large. The correlation between various cuts must be studied 
very carefully. 

In this section, we have only discussed the effect of lepton energy cuts 
on the observed ma distribution, and its dependence on unknown MSSM 
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parameters. It is possible that other cuts hke lepton isolation cuts, Efx cuts, 
and total Et cuts correlate with each other in a complicated manner to 
introduce additional uncertainties. This may be corrected quite easily by 
existing MC simulations. The aim of this section has been to point out an 
obvious source of uncertainty that has not been taken into account in current 
MC simulations, and to propose several observables that might constrain this 
uncertainty directly. 

4 Sensitivity to the Underlying Parameters 

In this section, we discuss if it is possible to extract the values of underlying 
MSSM parameters by measuring the mu distribution. As we have already 
seen the distribution depends strongly on m^-, and also on the X2Xi-^ and 
y^ll couplings. The latter is determined by the Higgsino components of X2 
and Xi- Because we take these neutralinos to be gaugino-like, their Higgsino 
components come from gaugino- Higgsino mixing. It depends on tan/5 and 
the Higgsino mass parameter /i, and |/i| is roughly equal to rriyXi, m^o or "111^+. 
Therefore the mu distribution gives at least one constraint on /x, tan/? and 
mj in addition to the well known constraint on m^o — m^o . 

I A,2 Ai 

Because there is no MC simulation including the parameter dependence 
of three body decay distributions, we estimate the sensitivity to tan 13 and 
mi under the following working assumptions: 

1. Backgrounds can be subtracted, or are negligible, and do not cause 
additional systematic errors. 

2. The dependence of the acceptance on mu can be corrected. 

Under these assumptions, we define the sensitivity function S as follows: 
S (Mi, M2, /i, tan/3, "i^-^^lfit; Mi, M2, /x, tan/5, mi^J^-,^^^^ 

= jE(^p-^r")'/^r"- (9) 

Here nf* {ui^^^ ) is the number of events in the i-th. bin of the mu distribution 
for the MSSM parameters (Mi, M2, /x, tan/3, m^- )|fit(input)- We normalize 

Y^i n^^ = Hi ^r^"* to some number N . Then S gives the deviation of the 
input distribution n^^^ from the distribution for the fit (nf*) in units of 
standard deviations. We takeQ N = 2500 and an mu bin size of 2 GeV. 

® The available number of lepton pairs from g decay ranges from 10^ to 1000 for 10 
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Figure 7: (a) Contours of constant sensitivity function S = 1, 2, 3, 4 in the (tan /3 
plane. See text for the definition of S . Input parameters are set (A) with mj = 250 GeV. 
For solid lines, we fix m^o, m^o and w~+ equal to those for parameter set (A) by varying 

Mf, Aff* and /i^', while tan/?*^* and mf^ are varied to see the sensitivity of the X2 decay 

distribution to these parameters. For the dot-dashed (dashed) lines, m%. = m'^^^ +(— )30 

X2 X2 

GeV have been taken. Diamonds correspond to the best fit points of each fit. (b) Contours 
of constant Br{x2 -^ Xi^~^^~) in the tan/3 and mj plane. The masses of X2; Xi; ^^id xt 
are set to equal to those for parameter set (A). 



In Fig. 7(a), we show contours of constant S = 1,2,3,4 (corresponding 
to la, 2a, 3a, Aa) in the (tan/^*^*, mp) plane. For the sohd hues, we take 
parameter set (A) and rrij = 250 GeV as input parameters, while for fitting 
parameters we vary tan/3 and mj = mj , fixing (Mf*, Mf*, /i^*) to reproduce 
the input values of (m^o, m^o, m-+).0 The resulting contours (solid lines) 
correspond to the sensitivity of the mu distribution to mj and tan/? when 
the three ino masses are known. 

In the figure, a strong upper bound on the slepton masses emerges, while 
the lower bound is weak. Notice that there is another minimum near tan f3 ~ 
2 and mj =170 GeV. By looking at the Af distribution, the two different 
minima may be distinguished. In that case a stronger lower bound would be 
obtained. On the other hand, mj < 260 (270) GeV is obtained if 5 < 1 (2) 
is required. This can be understood as a result of the large change of the 



fb^^ in the LHC study g] for a grand unified theory (GUT) scale gaugino mass parameter 
M = 100 GeV. At a high luminosity e'^e~ collider, 0(1000) lepton pairs could be produced. 
Our choice N — 2500 should thus not be too unrealistic. 

^In our fit, we fix the ino masses, therefore we implicitly assume that the "fake end 
point" problem discussed in section 2 can be resolved when the whole distributions are 
taken into account. 
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distribution between mj = 270 GeV and mj = 500 GeV found in Fig. 1. 
The mil distribution also constrains tan/? mildly. The constraint is not very 
strong due to our choice of parameters |/i| ^ M2; gaugino-Higgsino mixing 
is suppressed in this case. 

We computed these contours by fixing the three ino masses. This may 
not be a realistic assumption, because X2 might not be observed at the LHC. 
If we perform the fitting with varied m^\ , we can find almost the same mu 

'X.2 

distribution for different sets of parameters. In Fig. 7, the dashed (dot- 
dashed) lines correspond to constant S contours when m^}+ = m™+" — 30 

X2 X2 

GeV (m'!+ = rn^+^^ + 30 GeV) is required for the fitting. The best fit points 

X2 X2 

are tan/5 = 3.2 and mj = 233.6 GeV for m% = 772"+" — 30 GeV, and 

X2 X2 

tan/3 = 5 and m|- = 267 GeV for m^\ = ■02"+"* + 30 GeV. The best fit point 

X2 X2 

therefore moves from smaller tan j3 and mj to larger ones as m^\ is increased. 

X2 

Other distributions such as pip and Af might differ at different best fit points; 
however, we do not go into the details of these distributions. 

At each best fit point with different m%, 5 ~ 0. Therefore the dashed 
or dot-dashed contours are very similar to the contours of constant S where 
the parameters for the best fit point are taken as the inputs. It can be seen 

that as m^+ is decreased, the constraint on tan/3 becomes stronger, because 

X2 

Higgsino-gaugino mixing becomes more important. On the other hand, when 

m?+ is increased, the constraint on tan f3 tends to disappear. 

X2 

The constraint from the mu distribution is independent of that from the 
decay branching ratio. In Fig. 7(b), we show the contours of constant 
Br[x2 — > Xil~^l~). Notice that the branching ratio does not bound tan/3 
by itself unless Br{x2 ~^ Xi^~^^~) > 10%. The shapes of these contours are 
not particularly correlated to the constraint from the m,ii distribution shown 
in Fig. 7(a); therefore combining the two pieces of information might con- 
strain the parameter space further. However, we should be aware that the 
branching ratio depends on the whole sfermion mass spectrum. For example, 
if there is a substantial reduction of masses for third generation sfermions, 
the branching ratio will be changed. Fits using the m^u distribution are very 
important in that sense, because the shape depends only on ino mass param- 
eters and m,j. Notice also that we only measure the sum of the products of 
several branching ratios rather than Br[x2 ~^ Xi^^^) itself at the LHC. 
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5 Discussion and Conclusion 

In this paper, we examined the impact of neutrahno decay distributions 
on the study of the minimal supersymmetric standard model at the CERN 
LHC. The leptonic three body decay of the second lightest neutralino, X2 ~^ 
l~^l~Xii is known to be very important because the end point of the lepton 
invariant mass distribution m™^^ gives us direct information about the mass 
difference between X2 s-'^d Xij which gives us a stringent constraint on MSSM 
parameters. 

We found that the neutralino decay distribution depends on the slepton 
masses rather sensitively. Measuring the shape of the lepton invariant mass 
distribution can be important even for the determination of mfl'^^. In some 
cases, the measured end point may not coincide with the neutralino mass 
difference m^o — m^o, due to the strong suppression of the amplitude near 
m™'^^, while the leptonic branching ratio is still around a few %. 

On the other hand, if systematic errors can be removed, the distribution 
gives us a model independent probe of the slepton mass scale. For the X2 
signal from g decay, the S/N ratio is very large, or the background can be 
subtracted by using lepton pairs with opposite charge and different flavors. 
The remaining distribution can be further studied by looking at (A) the (3^o 
distribution measured by using events which have mu near its end point, and 
(B) the lepton transverse momentum asymmetry Af distribution, which is 
well correlated with the lepton energy asymmetry in the X2 '^^st frame. 

Notice that if the X2 momentum distribution is precisely measured, lep- 
tonic decay distributions may be discussed without any QCD uncertainty, 
which might otherwise be substantial. Of course, the correlation with other 
cuts (e.g., on I^t, J2 Et, lepton isolation) must be either small or determined 
from direct measurements, so that measurements in the lab frame allow us to 
reconstruct the X2 decay distribution in its rest frame. The revenue of such 
an effort to reduce the uncertainty from the cuts is information on slepton 
masses and neutralino mixings, independent of any assumption about the 
mechanism to break supersymmetry. 

In order to see if a study in such a direction is possible, dedicated MC 
simulations are necessary. Notice that the commonly available MC codes 
ISAJET and SPYTHIA used to simulate the three body decay distribution in 
phase space approximation. Our results show that a more careful treatment 
of X2 ~^ ^^^ Xi decays is important, if one is interested in decay distributions. 
Moreover, the acceptance of such di-lepton events depends on the lepton 
energy in X2 ""^ ^^^Xi decays. 
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The study of X2 decays can be easily done at e+e^ colliders. X2X2 P^-i^ Pro- 
duction or X2X1 co-production do not suffer from large backgrounds. Though 
the statistics is rather limited there, it would give us constraints on slepton 
masses and neutralino mixings. Notice that in supergravity models the lighter 
chargino and neutralinos are expected to be lighter than squarks and gluinos. 
Both LC and LHC may find and study X25 ^^^ the information gleaned from 
these analyses can be combined to obtain a better understanding of MSSM 
parameters and the SUSY breaking mechanism. 
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Appendix 

We show the explicit forms of the neutralino couplings used in section 2. 

The neutralino mass matrix Mn in the gauge eigenbasis {B, W^, H^, H^) 
is written as 



I Ml -mzcpsw mzSfsSw \ 

j^ ^ M2 mzc^cw -mzspcw 

-mzCfsSw mzCfsCw -/x 

V mzSfsSw -mzspcw -/i / 



(Al) 



where sw = sm9w, cw = cos9w, sp = sin/3, and C/j = cos/?. The mass 
eigenstates x° (« = 1 — 4) are related to the gauge eigenstates via the mixing 



matrix Nij as 



X° = NaB + N,2W^ + Ar,3^o + N,,Hl (A2) 



We neglect the possibility of CP violation for sparticles and take Nij as a 
real matrix by allowing negative values of my,o in Eq. (^ . 

The interaction Lagrangian of fermion-sfermion-neutralino and neutralino- 
Z^ couplings is written as 

Cint = -92Wa{xPL + b{xPR)fPx + (H.c.) + fzfMYl.x'AZl (A3) 
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where / stands for /, z/, d, and u. The couphngs (a, b) for leptons are, in the 
(Jl, Ir) basis, 

1 I TTl} I 

«AX = ^^ i [-^vl2 - NAltw]5x,L + Na3Sx,R \ , 

b'^X = -^\^^^NA35x,L + 2NAltw5x,R\, 

?^Ax = 0, (A4) 



where tw = tan6'vi/. We have always ignored 0{mi) terms in Eq. (^) in our 
numerical calculations. 

The coupling of x° and Z^ takes the form 

^BA = Nb^Naz - Nb4Na4. (A5) 
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